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ACTIVE  REMOTE  SENSING  USING  ULTRA  BROADBAND  RADIATION 


I.  Introduction 

The  possibility  of  detecting  the  constituents  of  the  atmosphere  and  measuring  atmospheric 
properties  was  proposed  soon  after  the  invention  of  the  laser.  LIDAR  (Light  Detection  and 
Ranging)  technology  encompasses  a  number  of  techniques  for  remote  sensing  of  the  atmosphere. 
The  fundamental  physical  processes  involved  in  lidar  include:  i)  Rayleigh  scattering,  ii)  Mie 
scattering,  iii)  Raman  scattering,  iv)  resonant  absorption,  v)  resonant  scattering,  and  vi) 
fluorescence.[l]  These  physical  processes  form  the  basis  of  lidar  techniques  that  can  be  used  in 
varied  environments.  Lidar  techniques  include:  i)  atmospheric  backscatter  lidar,  ii)  resonance 
scattering  lidar,  iii)  Differential  Absorption  Lidar  (DIAL),  iv)  fluorescence  lidar,  v)  Raman  lidar, 
and  vi)  Doppler  lidar. 

A  key  element  of  lidar  is  the  radiation  source.  Source  characteristics  tend  to  be  critical  in 
lidar  measurements  since  they  dictate  the  interaction  physics,  the  received  signal  and  ultimately 
the  signal  to  noise  ratio.  Relevant  source  characteristics  include:  i)  spectral  profile,  ii)  spectral 
brightness,  iii)  pulse  energy,  iv)  pulse  shape,  v)  pulse  duration,  vi)  pulse  repetition  rate,  vii)  beam 
divergence,  viii)  polarization,  and  ix)  coherence  properties. 

Active  remote  sensing  has  been  pursued  almost  exclusively  with  lasers  of  various  kinds. 
Lasers  are  highly  effective  in  many  respects  since  they  can  be  extremely  bright  and  highly  directed. 
In  many  applications,  however,  the  narrow  spectral  bandwidth  of  conventional  lasers  limits  their 
utility.  For  example,  in  DIAL  it  is  necessary  to  have  the  capability  to  scan  across  the  absorption 
spectrum  of  the  molecular  species.  In  many  cases  the  gain  bandwidth  of  the  available  source  is 
either  not  close  enough  to  the  molecular  transition  or  not  sufficiently  broad  to  encompass  the 
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absorption  spectrum.  It  is  desirable,  therefore,  to  have  a  source  that  shares  with  lasers  such  useful 
characteristics  as  brightness  and  directionality  and  yet  is  sufficiently  broadband  to  permit  target 
identification  over  a  wide  spectral  range. 

The  objective  of  this  report  is  to  discuss  the  remote  sensing  capability  of  ultra  broadband 
(UB)  sources  of  radiation  [2]  that  are  capable  of  detecting  and  monitoring  multi  species 
hazardous  air  pollutants  (HAPs).  Current  monitoring  techniques  are  not  capable  of  detecting 
many  of  the  HAPs  in  real  time  or  on  a  continuous  basis.  In  conventional  remote  sensing,  the 
received  signal  from  the  pollutant  is  usually  from  the  scattering  of  sunlight  or  a  single,  tunable 
laser  line.  The  relatively  narrow  tuning  ranges  of  lasers  and  the  inability  of  these  lasers  for  fast, 
real-time  tuning  restrict  the  use  of  LIDAR  and  DIAL  to  the  detection  and  ranging  of  a  select 
group  of  predetermined  pollutants. 

A  schematic  of  a  UB  remote  sensing  configuration  is  shown  in  Fig.  1.  A  modulated  laser 
interacting  with  a  nonlinear  medium  generates  UB  radiation  which  is  reflected  from  a  parabolic 
mirror  and  propagated  over  extended  distances  to  the  region  of  the  HAPs.  The  scattered 
radiation  having  the  spectral  signatures  of  the  HAPs  is  received  and  analyzed.  The  unique 
features  of  this  approach  are  the  high  average  power  possible  with  the  laser  beat  wave  UB  source 
and  the  capability  of  simultaneously  interrogating  multiple  HAPs. 

The  advent  of  ultrashort  laser  technology  has  ushered  a  new  era  of  time-resolved 
spectroscopy,  with  pulses  so  short  that  one  can  now  study  nonequilibrium  states  of  matter,  and 
explore  new  frontiers  in  science  and  technology.  One  of  the  most  important  nonlinear  optical 
processes  with  ultrashort  pulses  is  supercontinuum  generation  -  the  production  of  intense 
ultrashort  broadband  “white  light”  pulses.[2,3] 
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The  first  study  of  the  mechanism  and  generation  of  ultrashort  supercontinuum  dates  back 
to  1969,  when  Alfano  and  Shapiro  observed  the  first  “white”  picosecond  pulse  continuum  in 
liquids  and  solids.[2,4]  The  observed  spectra  extended  over  -  6000  cm*‘  in  the  visible  and 
infrared  wavelength  region.  The  large  spectral  broadening  of  ultrashort  pulses  was  attributed  to 
self-phase  modulation  (SPM)  arising  from  anharmonic  electronic  oscillations  and  a  model  of 
supercontinuum  generation  was  formulated.  Over  the  years,  advances  in  mode-locked  lasers  have 
led  to  the  production  of  wider  super-continua  in  the  visible,  ultraviolet,  and  infrared  wavelength 
regions  using  various  materials. 

The  supercontinuum  arises  from  the  propagation  of  intense  picosecond  or  shorter  laser 
pulses  through  condensed  or  gaseous  media.  When  an  intense  laser  pulse  propagates  through  a 
medium,  it  changes  the  refractive  index,  which  in  turn  changes  the  phase,  amplitude,  and 
frequency  of  the  pulse.  Supercontinuum  laser  sources  can  be  wavelength  selected  and  coded 
simultaneously  over  wide  spectral  ranges  (up  to  10,000  cm  *)  in  the  ultraviolet,  visible,  and 
infrared  regions  at  high  repetition  rates,  gigawatt  output  peak  powers,  and  femtosecond  pulse 
durations.  Ultrashort  supercontinuum  pulses  have  been  used  for  time-resolved  absorption 
spectroscopy  and  material  characterization.  The  applications  of  interest  here  include 
spectroscopy,  three-dimensional  imaging,  and  ranging  in  atmospheric  remote  sensing. 

The  plan  of  this  report  is  as  follows.  In  Section  n  we  discuss  the  generation  of  UB 
radiation.  Section  HI  presents  results  of  numerical  studies  of  UB  radiation  generation.  In  Section 
IV  we  show  that  with  the  use  of  three  wavelengths  a  significant  improvement  in  the  usual  (two- 
wavelength)  DIAL  method  is  possible.  In  particular,  it  is  shown  that  knowledge  of  the 
dependence  of  the  reflectivity  and  extinction  coefficient  on  wavelength  is  not  necessary  if  three  or 
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more  wavelengths  are  used.  Section  V  presents  a  discussion  of  an  experiment  where  UB 
radiation  is  generated  in  the  near  IR,  with  application  to  the  detection  of  HAPs.  Section  VI 
summarizes  the  report. 

n.  Ultra  Broadband  Radiation  Generation 

The  purpose  of  this  section  is  to  present  an  outline  of  the  analysis  that  demonstrates  UB 
radiation  generation  from  a  single  pulse  or  from  the  beating  of  two  laser  pulses  with  differing 
frequencies  (laser  beat  waves).  Ultra  broadband  radiation  can  be  generated  in  a  nonlinear  optical 
medium  with  pico-second  laser  pulses.  In  the  past,  the  single  pulse  conversion  efficiency  and 
bandwidth  have  been  impressive,  but  average  power  has  been  limited  to  relatively  modest 
values,[2,3]  Recently,  a  new  method  has  been  proposed  for  generating  high  average  power  UB 
radiation.  [5]  In  this  method,  UB  radiation  is  generated  by  beating  two  laser  beams  with  slightly 
different  frequencies  in  a  nonlinear  medium.  [5]  The  conversion  efficiency  for  this  process  can  be 
shown  to  be  relatively  high,  and  the  bandwidth  of  the  radiation  can  be  on  the  order  of  100%. 
Especially  in  the  3-8  ptm  wavelength  region,  where  there  is  a  lack  of  available  laser  sources,  the 
UB  radiation  can  provide  the  necessary  illumination  required  for  active  remote  sensing.  The 
source  size  of  the  UB  radiation  is  extremely  small,  which  allows  for  beaming  the  radiation 
extended  distances,  beyond  3  km. 

Lasers  with  optical  (-  0.5  |xm)  and  near-IR  (~  1  p-m)  wavelengths  can  be  used  to  generate 
UB  from  -  0.4  pm  to  ~  3  pm.  To  reach  longer  wavelength  regions,  CO2  lasers  may  be  used. 
These  lasers  are  line-tunable  in  the  9-1 1  pm  range.  If  two  closely-spaced  CO2  laser  lines  are 
selected,  the  beating  of  these  lines  in  an  appropriate  nonlinear  medium  such  as  gallium  arsenide 
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can  lead  to  substantial  spectral  broadening,  thus  covering  the  mid  and  far  IR  regions  of  the 
spectrum.  In  principle,  the  conversion  efficiency  and  bandwidth  can  be  optimized  by  selecting  the 
appropriate  nonlinear  medium  and  adjusting  the  beat  frequency,  power,  and  spot  size  of  the  laser. 

A.  UB  Generation  Using  a  Short  Laser  Pulse 

The  longitudinal  gradients  in  the  intensity  of  a  single  laser  pulse  can  lead  to  frequency  up- 
shifting  and  broadening  via  self-phase  modulation.  Consider  the  nonlinear  index  of  refraction 
lesulting  from  the  (third-order  susceptibility)  nonlinearity,  n  =  no  +  n2l,  where  no  is  the  linear 
index,  I  is  the  laser  intensity,  n2  is  the  nonlinear  index,  n2l  =  27CX^^EVno,  and  E  is  the  electric  field 
of  the  laser  pulse.  The  change  in  phase  5<|)  for  a  laser  field  E  ~  expi(koz  -  tOjt  +  6<|))  which 
propagates  into  a  nonlinear  medium  is  given  by  8<|)  =  (tOo/c)  J  dzn2l  =  (Gb/c)n2l2,  where  (Oj  and  ko 
are  the  nominal  frequency  and  wave  number  of  the  laser  field  and  z  is  the  axial  propagation 
distance.  The  change  in  frequency  is  given  by  5©  =  -35<|)/3t.  Hence,  as  the  laser  pulse  propagates, 
new  frequencies  are  generated  [2,6,7] 

5©  =  -z(©o  /  c)n29l  /  9t.  (1) 

In  particular,  at  the  front  of  the  pulse  dVdt  >  0  and  the  local  frequency  is  down-shifted,  whereas  at 
the  back  of  the  pulse  81/3 1  <  0  and  the  local  frequency  is  up-shifted.  The  above  simple  estimate  is 
valid  for  small  frequency  shifts  |8a)| «  ©q.  For  ultra-short  pulses,  the  intensity  gradients  can  be 

very  large  which  results  in  large  frequency  shifts  and  the  generation  of  a  supercontinuum. 

Spectral  broadening  has  been  observed  experimentally  with  ultrashort  laser  pulses,  i.e., 
pulse  durations  on  the  order  t  <  1  ps.  Using  femtosecond  laser  pulses  with  a  wavelength  Xo  ~  1 
|xm,  a  spectral  broadening  of  20,000  cm  '  is  readily  achievable  in  a  condensed  medium.[2]  Short 
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pulses  are  necessary  in  order  to  achieve  large  intensity  gradients,  3l/3t  ~  I/t.  The  use  of  ultrashort 
pulses  has  its  limitations:  (i)  the  supercontinuum  generated  is  also  short  pulsed  (<  1  ps)  and, 
perhaps  more  importantly,  (ii)  the  repetition  rate  and  hence  the  average  power  of  the  generated 
radiation  is  low,  due  to  the  low  repetition  rate  and  low  average  power  of  the  short  pulse  laser. 


B.  UB  Generation  Using  Two  Beating  Lasers 

Supercontinuum  generation,  in  principle,  can  be  achieved  using  the  beat  wave  produced  by 
the  interaction  of  two  long  pulse  laser  beams  of  two  different  frequencies,  (Oi  and  (p2,  within  a 
nonlinear  medium.  Consider,  for  example,  two  Unearly  polarized  lasers  with  electric  fields  given 

by 

El  =  El  cos(kiZ-(Oit)ex, 

E2  =  E2  cos(k2Z  -  0)2 


where  ki.2  are  the  wave  numbers  and  Cj  is  the  unit  vector.  Within  the  nonlinear  medium,  these 
two  fields  will  interact  to  produce  a  beat  wave  via  the  term  in  the  nonlinear  polarization 
field.  Assuming  equal  amplitudes,  Ei  =  E2  =  Aq  /  2,  we  have 


E  =  Ei+ E2  =  Aq  cosfkoz  - 


cDqOcos 


(f- 


(3) 


where  0)1,2  =  teo  ±  Ao)/2,  ki.2  =  ko  ±  Ak/2,  and  0)i  -  ojj  =  Am.  Hence,  |Ei  +  E2P  represents  a  beat 
wave,  i.e.,  the  intensity  envelope  I  of  the  two  laser  fields  is  modulated  at  the  beat  frequency  Ao) = 
0)|-0)2, 
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I  ~  ^|Ei  +  E2P ^  =  (Aq  /  4)[1  +  cos(Akz -  Acot)], 


(4) 


where  <  >  represents  a  time  average  over  the  fast  frequency  (Oo-  In  effect,  the  longitudinal 
intensity  profile  resembles  a  train  of  short  pulses  of  duration  t  ~  1/Aco  and  results  in  large 
gradients  in  the  intensity  3I/3t  ~  Ami.  Likewise,  the  induced  frequency,  as  given  by  Eq.  (1), 

8ci)  /  (Dq  -  z(Aco  /  c)n2l,  (5) 

can  be  large.  As  an  example,  consider  two  laser  beams  operating  on  two  CO2  lines,  =  10.5  |im 
and  =  10.6  ^im.  Hence,  Aco  =  litcilfki-H'ki)  =  1.7  x  10*^  sec'*  and  x  =  1/Am  =  590  fs.  Since 
the  two  laser  beams  are  long  pulse,  the  resulting  radiation  can  be  long  pulse  with  relatively  high 
average  power. 


C.  Fully  Nonlinear  and  Dispersive  Model  of  Spectral  Broadening 

The  two-dimensional  nonlinear  evolution  of  a  short  laser  pulse  in  a  nonlinear  medium  is 


described  by  the  wave  equation 


c^  dt^ 


J 


4n 


at 


Pn, 


2  -P 


(6) 


where  k j.  =  2/ro  is  the  transverse  wave  number,  ro  is  the  spot  size,  and  Pp  is  the  complete 


polarization  field.  It  is  convenient  to  use  the  notation  E  =  A  expj^i(koZ  —  mot)jex ,  where  A  is  the 


slowly  varying  electric  field  amplitude.  In  the  following,  kj^  =  0  will  be  assumed. 

Injecting  two  laser  beams  of  slightly  different  frequencies  into  the  nonlinear  medium  can 
provide  spectral  broadening  together  with  high  average  power.  The  initial  optical  field  is  taken  to 
be 
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E(z  =  0)  =  (Aq  /  2)cos((ko  -  Ak  /  2)z  -  (cOq  -  Aod  /  2)t)ex 


+  (  Aq  /  2)  cos((ko  +  Ak  /  2)z  -  (coq  +  Ato  /  2)t)ex 
=  Aq  cos(A(ot  /  2)  cos(a)o'^)ex ,  (7) 

where  Ao  is  the  real  and  constant  amplitude,  ko  =  nooWc  =  (Oo/vg,  Ak  =  Aco/Vg  is  the  difference  in 
wave  numbers  of  the  two  lasers,  Aco  is  the  frequency  difference,  x  =  t  -  z/vg,  Vg  is  the  linear  group 


velocity.  The  power  associated  with  the  optical  field  is  P  =  cro  no|Ap  / 16 .  It  is  convenient  to 


^1/i^  /  9  \l/2 

introduce  the  complex  field  amplitude  P  =  (cro  no  /  IhPjsj  1  A,  such  that 


pl/2 


=  P  /  Pn  is 


the  normalized  power,  where  P^  =  /  (o)otto*'2 )  critical  power  for  nonlinear  self- 

focusing.  Defining  z  =  z  /  Lq,  Lq  =  noZR  /  2,  and  Zr  =  (OoTo  /  2c,  the  wave  equation  can  be 
solved  for  the  complex  amplitude  and  cast  into  the  following  normalized  form: 
ap*^^  /  az = i(i  -I-  iicooTor^d  /  ax)(|pp'^) 

-  i(i  /  2)R2a^p’'^  /  ax^  -t-  (1  /  6)R3a^p‘'^  /  ax^  (8) 


Here  x  =  x  /  Xq  =  (t  -  z  /  Vg )  /  Xq  is  the  normalized  time  in  a  frame  moving  with  the  group 

velocity,  Xo  is  a  characteristic  pulse  time  or  period,  R2  =  Lq/Ldi  and  R3  =  WLds-  The  first  and 
second  order  dispersion  distances  are,  respectively, 

LD2=t§/H2,  (9a) 

L-DS  ~ '*’0  ^  ^3’ 
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where  in  the  kj.  =  0  limit,  the  quantities  H2  and  H3  are  the  usual  group  velocity  dispersion  (GVD) 
coefficients  P2  and  P3  [3].  For  P  =  Pn,  the  distance  Lo  is  the  same  as  the  “nonlinear  distance” 
described  in  Ref.  3.  The  initial  value  of  P^^^  at  z  =  0  is 

P''^(z  =  0,T)  =  (Po(0,T)/  Pn)’'^  cos(Acdt/  2) ,  (10) 

where  Po(0,x)  is  the  initial  axial  profile  of  the  power. 

In  the  absence  of  dispersion,  R2  =  R3  =  0,  the  reduced  wave  equation  in  Eq.  (8)  becomes 


*-»*  1  ^  1/2 

The  evolution  of  the  amplitude  Pq  and  phase  (j)  of  P  is  given  by 

—  -P„  ,  (12a) 

||=(i-(os‘a<p/at)Po''^  ci2b) 

where  P’^^  =  Pq^^  exp(i<l>)  and  the  initial  conditions  (at  z  =  0)  for  Pq^^  are  given  by  Eq.  (10). 

Equations  (8)  can  be  used  to  describe  self-phase  modulation,  continuum  generation,  and 
soliton  formation,  including  the  effects  of  group  velocity  dispersion  and  finite  The  standard 
theory  [2]  of  self-phase  modulation  and  continuum  generation  for  ultrashort  pulses  neglects  the 
effects  of  group  velocity  dispersion,  i.e.,  the  terms  proportional  to  R2  and  R3  are  neglected.  In 
standard  soliton  theory  [3,8],  the  first  term  on  the  right  of  Eq.  (8)  is  neglected.  Equation  (8) 
allows  continuum  generation  and  solitons  to  be  studied  self-consistently  including  these  effects. 
In  particular,  it  can  be  used  to  study  self-phase  modulation,  continuum  generation,  and  soliton 
formation  for  a  beat  wave  produced  by  two  lasers  of  different  frequencies,  as  described  by  Eqs. 
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(2a)  and  (2b).  Numerical  solutions  to  Eq.  (8)  are  described  in  the  next  section.  These  results 
show  that  dispersion  effects  can  significantly  modify  the  spectrum  even  when  the  dispersive  terms 
R2  and  R3  are  small. 


D.  Lowest  Order  Solution  in  the  Absence  of  Dispersion 

Approximate  solutions  to  Eq.  (12),  for  small  variations  in  the  amplitude  and  frequency,  are 


P  =  PoCOS^(A(DT/2) 


2  (Oq  Pn  Lo  } 


(13a) 


©  =  (Oo 


,  1  Ao)  Pq  z  .  . .  . 

1+ - - sin(A(DX) 

2  coo  Pn  Lq 


(13b) 


where  A©  =  ojdAX/Xo  is  the  initial  frequency  difference,  and  AX  is  the  wavelength  difference. 


III.  Numerical  Results  for  UB  Radiation  Generation 

Equation  (8)  has  been  solved  numerically  using  the  split-step  Fourier  method,  a 
pseudospectral  numerical  technique  which  has  been  widely-used  in  modeling  nonlinear  effects  in 
fiber  optics.[3,9]  The  initial  amplitude  P^'^  (z  =  0,  x)  is  specified,  and  the  complex  ampUtude 
P^^^  is  calculated  using  a  step  size  Az  which  is  small  compared  with  Lo-  Since  a  fast  Fourier 
transform  (FFT)  is  used,  the  number  of  grid  points  in  t  must  be  a  power  of  two,  and  the 
waveform  must  go  to  zero  near  the  boundaries  to  avoid  aliasing  effects.  For  the  two  laser  beat- 
wave  case,  it  is  convenient  to  chose  Xo  =  A©”'  and  specify  Po  =  Pn.  In  addition,  one  must  specify 
the  normalized  input  parameters  A©/©o,  Ld2/Lo,  and  Loa/Lo. 
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Ultra  broadband  radiation  is  traditionally  produced  by  using  a  single  short  pulse  or  a  train 
of  isolated  pulses,  and  numerical  solutions  to  an  equivalent  form  of  Eq.  (8)  for  such  pulses  are 

extensively  described  in  Ref.  3.  Figure  2  plots  the  normalized  power  P  (x)  at  locations  z/Lo  =  0 
and  17  for  a  single  pulse  with  (z  =  0,  x)  =  exp(-x^/2Xo^)  and  tObXo  =  45.  This  Gaussian  pulse 
case  neglects  dispersion.  One  can  show  analytically  [3]  that  in  the  limit  of  zero  dispersion,  the 
steepening  on  the  back  side  of  the  pulse  should  reach  an  infinite  slope  at  z/Lo  =  0.39cabXo,  a  result 
which  is  consistent  with  the  numerical  solution  in  Fig.  2.  This  steepening  causes  substantial 
broadening  in  the  power  spectrum  S(©),  which  is  calculated  by  taking  S(a))  =  P|J^^(0))  ( P]J^^(©) ) 
where  is  the  Fourier  transform  of  P^^^(x),  and  (Pk^^(a)))*  is  its  complex  conjugate. 

Figure  3  plots  S(©)  at  z/Lo  =  17  and  shows  the  frequency  spread  is  comparable  to  Cflb,  indicating 
the  generation  of  ultra  broadband  radiation.  The  oscillations  in  the  spectrum  are  consistent  with 
the  modulation  commonly  observed  in  many  experiments.[2,10] 

Qualitatively  similar  spectra  can  be  produced  for  the  laser  beat- wave  technique  described 
in  Section  11.  Figure  4  shows  the  power  spectrum  S((i))  at  z/Lo  =  80  for  a  case  with  A(o/C0o=  0.01 
and  a  very  small  amount  of  dispersion  (R2  =  3x10"^).  The  power  spectrum  consists  of  discrete 
lines  separated  by  Am  and  shows  signiHcant  spectral  broadening,  especially  in  the  anti-Stokes 
direction  (©  >  coo).  Even  this  small  amount  of  dispersion  is  sufficient  to  change  the  numerical 
solution  sigiuficantly.  For  Rj  =  0,  the  model  gives  an  optical  shock  and  substantial  spectral 
broadening  at  z/Lo  =  45,  and  the  numerical  method  breaks  down  shortly  after  that  point. 

Although  most  analyses  of  ultra  broadband  radiation  neglect  dispersion,  it  is  clear  that 
dispersion  plays  a  major  role  when  parameters  corresponding  to  most  real  materials  and 
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experiments  are  used.  For  a  single  Gaussian  pulse  with  a  small  but  nonzero  R2,  the  pulse  flattens 
and  spreads  after  propagating  in  z  over  a  distance  which  is  typically  only  10-20%  of  the 
characteristic  dispersion  distance  Ld2.[3]  This  is  accompanied  by  a  strong  steepening  at  the  edge 
of  the  pulse  and  the  generation  of  high  frequency  edge  oscillations.[3]  An  analogous  process  is 
predicted  to  occur  in  the  beating  laser  case,  except  that  the  spreading  of  the  waveform  on  the  Aco 
time  scale  is  clearly  constrained  by  the  next  half-oscillation  in  the  waveform. 

Figure  5  plots  P  (x)  at  z/U  =  0, 10, 20,  and  30  for  a  beating  laser  case  with  Aco/tOo  =  0.01, 
R2  =  3.8x10 and  R3  =  4x10"*.  These  parameters  are  appropriate  for  the  proof-of-principle 
experiment  described  in  Section  V,  operating  at  Pq  =  1  in  a  fused  silica  optical  fiber  with  a  peak 

intensity  of  2xl0"  W/cml  The  figure  shows  three  half-oscillations  of  the  beating  waveform.  As 
the  laser  propagates,  fast  oscillations  eventually  form  at  the  boundaries  between  successive  half¬ 
oscillations.  Figure  6  plots  the  power  spectrum  S(co)  at  z/Lo  =  30.  TTie  spectrum  is  relatively  flat 
except  for  the  harmonic  structure  at  Am  and  extends  to  about  15%  of  the  central  frequency  CH)  in 
both  the  Stokes  and  anti-Stokes  directions.  This  spread  is  typical  of  many  single  pulse 
experiments  in  optical  fibers.[2,3,10] 

Larger  spreads  in  the  power  spectrum  can  in  principle  be  produced  by  optimizing  the 
frequency  mismatch  Am,  increasing  the  laser  intensity,  propagating  further  into  the  fiber,  or 
choosing  a  material  with  a  lower  dispersion  parameter  at  1.06  flm.  Higher  laser  intensities  can 
lead  to  breakdown  and  permanent  damage  to  the  material,  while  propagating  too  far  into  the 
material  produces  little  additional  spreading  and  tends  to  give  a  nosier  spectrum.  The  choice  of 
the  optimal  Am  involves  a  tradeoff  between  the  propagation  distance  required  for  the  desired 
frequency  spread  and  dispersion  effects.  Equation  (13a)  shows  that  reducing  Aoi/mo  can 
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substantially  reduce  undesirable  dispersion  effects,  but  Eq.  (13b)  shows  that  it  increases  the 
distance  z  required  to  achieve  a  given  spread  in  oycoo- 

The  numerical  model  has  also  been  used  to  evaluate  potential  ultra  broadband  sources 
using  two  closely-spaced  lines  from  a  CO2  laser  and  gallium  arsenide  as  the  nonlinear  medium. 
Again  dispersion  plays  a  major  role  in  determining  the  spectrum.  Gallium  arsenide  exhibits 
anomalous  dispersion  (P2  <  0)  in  this  regime,  which  results  in  cusp-like  spikes  on  the  waveform 
and  significant  spectral  broadening  after  only  a  few  Rayleigh  lengths  of  propagation.  Preliminary 
results  suggest  that  in  the  9  |j.m  regime,  a  2%  mismatch  in  the  wavelengths  may  produce 
significant  radiation  down  to  6  |im.  This  case,  which  is  shown  in  Fig.  7,  used  Lo/Ld2  =  0.2  and 
L(/Ld3  =  0.027,  which  corresponds  to  P  =  Pn  tmd  an  intensity  of  10*°  W/cm^.  Spectral  broadening 
can  be  overwhelmed  by  dispersion  when  the  dispersion  length  is  comparable  to  the  Rayleigh 
length,  so  relatively  small  changes  in  parameters  could  modify  the  spectrum  considerably  in  this 
parameter  regime.  Ultra  broadband  radiation  using  a  single  pulse  CO2  laser  was  first  reported  by 
Corkum,  et  al.[l  1] 

rv.  Three  Wavelength  DIAL 

The  use  of  DIAL  for  remote  sensing  of  chemical  species  in  the  atmosphere  is  a  well 
known  technique  that  can  detect  and  identify  a  wide  range  of  molecular  constituents.  However, 
DIAL  is  presently  operated  with  two  distinct  laser  wavelengths  when  the  absorption  coefficients 
of  the  chemical  species  are  known  as  a  function  of  wavelength.  When  using  only  two  frequencies, 
one  needs  to  know  the  target  reflectivity  and  the  extinction  coefficient  of  the  atmosphere  (which 
are  not  typically  known)  at  these  wavelengths  to  obtain  accurate  measurements. 
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Knowledge  of  the  dependence  of  the  reflectivity  and  extinction  coefficient  on  wavelength 
is  not  necessary  if  three  or  more  frequencies  are  used.  If  a  third  frequency  is  used,  which  is 
symmetrically  located  on  either  side  of  the  absorption  peak,  as  shown  in  Fig.  8,  it  can  be  shown 
that  the  pollutant  density  can  be  determined  independently  of  the  reflectivity  and  extinction 
coefficient  of  the  background.  We  are  assuming  here  that  the  dependence  of  the  background 
reflectivity  and  extinction  coefficient  have  a  weak  linear  dependence  on  wavelength.  The 
pollutant  density  is  then  simply  a  function  of  the  received  power  at  the  three  wavelengths.  Using 
three  specific  wavelengths  can  therefore  yield  far  mote  accurate  values  for  the  pollutant  density. 
The  UB  transmitting  signal  can  serve  as  a  three  wavelength  DIAL  system  over  multiple 
absorption  bands.  Each  pollutant  will  absorb  its  characteristic  radiation,  which  can  be  detected, 
identified,  and  analyzed,  yielding  both  the  range  and  concentration. 

Using  the  LIDAR  equation  [1]  it  can  be  shown  that  using  three  wavelengths,  Xo,  Xi,  and 

X2,  the  minimum  density  of  the  pollutant  N  is  given  by 

<n(VP,(>.i)P,(>.2)/P,(>.o)) 

2R(o(>.o)-o(\i)) 

where  Pr(X)  is  the  power  received  at  the  detector  at  wavelength  X,  R  is  the  range,  and  a(X) 
denotes  the  absorption  cross  section  at  wavelength  X.  As  an  illustration  consider  a  pollutant  at  a 
range  of  R  ~  1  km  and  a  signal  strength  Pr(Xi,2)/Pr(Xo)  “  2.  The  sensitivity  of  a  three  wavelength 
DIAL  system  is  N/No  =  1.3  x  lO"*  for  the  pollutant  SO2  having  a  cross  section  of  a(X  =  0.8-1 
|im) «  10*'’  cm^,  where  No  is  the  ambient  air  density. 

V.  Experimental  Approach  for  Ultra  Broadband  Radiation  Generation 
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The  generation  of  UB  radiation  will  be  carried  out  at  the  Naval  Research  Laboratory 
(NRL)  using  available  laser  facilities.  Single  laser  pulses  with  sub-picosecond  pulse  durations 
from  the  Table  Top  Terawatt  laser  have  been  used  to  demonstrate  the  generation  of  ultra 
broadband  radiation.  High  average  power  UB  radiation  in  nonlinear  media  can  be  generated  with 
long  pulse  lasers  modulated  at  subpicosecond  periods.  This  can  be  achieved  in  two  possible 
arrangements.  The  first  is  to  inject  two  laser  pulses  with  slightly  different  frequencies  in  a 
nonlinear  medium,  and  the  second  is  to  inject  a  single  laser  pulse  into  a  Raman  cell  in  order  to 
generate  a  second  frequency  within  the  pulse  prior  to  injection  into  the  nonlinear  medium. 

Two  mode-locked  lasers  at  NRL  are  available  for  this  experiment.  One  laser  is  a  Nd:YAG 
laser  with  wavelength  equal  to  1.064  |xm,  energy  of  300  mJ  per  pulse  and  pulse  duration  of  300 
psec.  Another  laser  is  a  Nd:Glass  laser  with  wavelength  equal  to  1.053  |im  and  similar  pulse 
energy  and  duration.  These  two  lasers  can  be  synchronized  with  each  other  by  using  the  same  RF 
driver  for  their  mode-lockers.  Exact  synchronization  will  be  achieved  with  optical  delay  lines. 

The  laser  beams  are  then  adjusted  to  have  the  same  divergence,  so  that,  when  focused  they  will 
overlap  at  the  same  axial  position.  The  lasers  will  also  be  aligned  to  overlap  transversely  at  the 
focal  plane.  The  angular  separation  of  the  two  lasers  will  also  be  minimized  to  obtain  the 
maximum  overlap.  The  wavelength  difference  of  10  nm  between  the  two  lasers  will  create  a 
modulation  period  of  300  fsec.  The  pulse  energy  of  the  two  lasers  will  be  adjusted  to  be 
approximately  the  same  to  ensure  large  amplitude  modulation  in  the  beat  wave  created. 

Alternatively,  a  single  laser  may  be  used  to  generate  the  modulated  pulse.  A  Raman  cell 
containing  a  suitable  nonlinear  medium  can  be  used  to  generate  a  second  wavelength  utilizing  the 
Raman  interaction.  Different  nonlinear  media  can  impart  different  Raman  shifts  from  the 
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fundamental  for  the  second  laser  frequency.  For  example,  carbon  disulfide  can  produce  a  Raman 
shift  of  ~80  nm  for  a  fundamental  wavelength  of  1.064  fim,  while  carbon  tetrachloride  will 
generate  a  shift  of  ~56  nm.  These  wavelength  differences  between  the  two  laser  frequencies  wUl 
generate  modulation  periods  of  38  fsec  and  54  fsec  respectively.  The  advantage  of  using  the 
Raman  cell  to  generate  the  second  laser  frequency  is  that  the  two  laser  lines  are  naturally 
synchronized  and  aligned.  However,  it  is  more  difficult  to  adjust  the  relative  energies  in  the  two 
laser  lines  by  controlling  the  conversion  efficiency.  One  could  split  the  two  laser  lines  and 
attenuate  one  of  them  to  obtain  equal  proportion  of  energy  in  each  line.  However,  the  necessity 
of  recombining  the  two  lasers  requires  delicate  alignment  similar  to  using  two  different  lasers. 

The  modulated  laser  beam  will  then  pass  through  a  suitable  nonlinear  medium  for 
generating  the  UB  radiation.  The  sub-picosecond  structures  in  the  laser  pulse  will  generate  the 
UB  radiation  through  the  mechanism  of  self-phase  modulation.  Experiments  will  be  performed 
using  various  media  with  high  nonlinear  index  of  refraction  including  liquids  and  solid  state 
materials.  Liquids  have  the  advantage  of  self-healing,  which  is  desirable  for  laser  intensities 
approaching  the  laser  breakdown  thresholds.  Solid  material  will  be  used  at  lower  laser  intensities. 
Interaction  distances  long  compared  to  the  diffraction  length  are  desirable,  and,  hence,  optical 
fibers  are  natural  candidates  for  the  nonlinear  medium.  Suitable  optical  fibers  will  be  used  to 
study  the  UB  radiation  generation  by  modulated  long  laser  pulses. 

VI.  Summary 

Ultra  broadband  (UB)  radiation  is  a  promising  source  for  active  remote  sensing.  It  can  be 
readily  produced  by  propagating  intense,  short  laser  pulses  in  appropriate  solids  or  liquids.  In  this 
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report,  we  assess  the  feasibility  of  using  two  long  pulse,  mode-locked  lasers  with  a  small 
separation  Aco  in  frequency  to  generate  UB  radiation.  Such  a  system  may  be  less  expensive  and 
more  compact  than  single  rep-rated  short  pulse  lasers  while  operating  at  substantially  higher 
average  powers.  The  theory  of  UB  generation  with  laser  beat  waves  is  presented,  including  the 
effects  of  self-phase  modulation,  nonlinear  steepening,  and  both  first  and  second  order  dispersion. 
Numerical  solutions  to  the  model  equation  show  parameters  regimes  where  UB  generation  may 
be  expected.  Dispersion  plays  a  major  role  in  the  spectral  broadening,  so  both  laser  and  material 
parameters  must  be  carefully  chosen.  As  expected,  the  behavior  for  laser  beat  waves  is  similar  in 
many  cases  to  that  for  a  single  pulse  with  the  pulse  width  replaced  by  the  beat  wave  period,  Aco  ^ 

One  application  of  the  UB  source  is  for  multi-wavelength  DIAL  to  determine 
concentrations  of  hazardous  pollutants.  Conventional  DIAL  systems  use  only  two  wave-lengths, 
but  the  addition  of  a  third  wavelength  can  eliminate  the  requirement  that  the  reflectivity  and 
extinction  coefficients  be  known. 

Finally,  a  proof-of-principle  experiment  for  UB  generation  using  laser  beat  waves  is 
outlined.  The  experiment  would  use  two  existing  mode-locked  lasers  at  NRL  operating  at  1.064 
Hm  and  1.053  jim,  with  an  alternative  approach  using  a  single  laser  and  a  Raman  cell  to  shift  a 
portion  of  the  radiation  for  the  second  wavelength.  Initial  experiments  will  be  carried  out  using 
liquids  or  optical  fibers  as  the  nonlinear  media. 
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Rg.  1.  Schematic  of  active  rmote  sensing  using  ultra  broadband  radiation- 
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1.5 


Hg.2.  Plot  of  the  nomalized  power  P  (t)  =  P(x)/Pn  at  2/Lo=0  (dotted  line)  and  17  (solid) 
for  a  single  pulse  with  a  Gaussian  shape  P^'^(z  =  0,  t)  =  exp(-t*/2Xo^)  and  a^Xo  = 
45,  This  case  neglects  dispersion  and  is  calculated  nuraaically  from  Eq.  (12). 
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Plot  of  the  power  spectrum  S((d)  in  arbitrary  units  at  z/Lo  -  17  for  the  Gaussian 
pulse  case  shown  in  Fig.  2,  The  frequency  spread  is  comparable  with  the  initial 


laser  frequency. 
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Rg.  4.  Model  calculation  of  the  power  spectrum  (arbitraiy  units)  for  a  laser  beat  wave 
example  with  frequency  separation  Aca/Oo — 0.01  and  the  dispersion  parameter  Rj 
0.0003.  Spectrum  is  similar  to  the  single  pulse  case  but  shows  discrete  lines 
separated  by  Am. 
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6.  Plot  of  the  calculated  power  spectrum  S(ei))  in  arbitrary  units  at  z/Lo  =  30  for  the 


laser  beat-wave  example  shown  in  Fig.  5. 


Rg.  7.  Plot  of  the  calculated  power  spectrum  S((o)  at  z/Lo = 4  for  a  case  with  Aaifah  = 

0.02,  R2  =  -0.20,  and  R3  =  0.027.  Tlie  parameters  for  this  anomalous  dispersion 
case  (R2  <  0)  are  appropriate  for  using  two  lines  from  a  CC)2  laser  as  the  source  and 
gflllinm  arsenide  as  the  nonlinear  medium 
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Hg.  8.  Locations  of  tiuee  wavelengths  around  an  absorption  peak.  The  differences  in 
received  power  are  used  in  a  three  wavelength  DIAL  scheme  to  calculate 
concentrations  of  the  trace  specie. 
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